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ABSTRACT 


The equations developed by Riehl (1963) for a steady, 
symmetrical hurricane are used to specify a vortex over a 
fixed depth boundary layer, which interacts with a simple 
ocean model. Only vertical mixing in the surface layer is 
included in the oceanic portion of the model. Interaction 
of the atmosphere and the ocean occurs through the predic- 
tion of boundary layer parameters. Bulk aerodynamic equa- 
tions are used to calculate the air-sea interaction. 
Comparison of the predicted boundary parameters with ob- 
served mature tropical storms shows good agreement for all 
variables except the specific humidity, which exceeas 
saturation. Two attempts to correct for this situation have 
shown that small changes in the surface moisture and tempera- 
ture values strongly affect the maximum wind speed predicted. 
The proportion of the water vapor transported out of the top 
of the boundary layer without adding latent heat to the 
boundary layer appears to be an important factor in deter- 
mining the intensity of a steady-state hurricane. 

Although the ocean model is quite simple, realistic 
horizontal oceanic temperature profiles with respect to the 
hurricane center were obtained. However the magnitude of 
changes in oceanic heat content and the sea-surface tempera- 
tures were smaller than expected in actual hurricanes, since 
advection and upwelling Saee were not included in the 


ocean model, 
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I, INTRODUCTION 


Observational and numerical investigations of tropical 
hurricanes have been made by many tropical meteorologists. 
Dynamical and structural features of mature three-—dimensional 
tropical storms are not completely clear but much has been 
learned through the use of steady-state axisymmetric hurri- 
cane models. Malkus and Riehl (1960) examined the dynamics 
and energy transformations in such a steady-state hurricane. 

Investigation of inflow trajectories, shown in Figure l, 
gave insight to the oceanic heat source necessary for the 
maintenance of a steady-state storm. The dynamic model of 
the inflow layer developed useful relations between wind 
speed, pressure gradient, surface shearing stress, mass flow 
and convergence. 

In a more recent work, Riehl (1963) examined some rela- 
tions between wind and the thermal structure of the steady- 
State hurricane. The steady-state hurricane is defined as 
one in which a balance exists between the heat source and 
momentum sink at the air-sea interface. In terms of a bulk 
aerodynamic estimate of these processes, the heat source de- 
pends linearly on wind speed while momentum sink depends on 
the square of the wind speed. A steady hurricane model was 
developed to compare with observed storms. A two-layer model 
was considered with conservation of absolute angular momentum 
assumed in the outflow layer, shown in Figure 1, and conser- 


vation of potential vorticity assumed in the inflow. With 
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these assumptions the momentum transport at the air-sea 
interface can be specified, and thus the radial profile of 
the tangential wind component. By assuming the kinetic ener- 
gy necessary to drive a hurricane is provided by the oceanic 
heat source, Riehl estimated the heat source required for a 
given wind field to exist in steady state. 

The purpose of this study is to utilize the steady-state 
hurricane model proposed by Riehl in experiments with an 
interacting ocean model. Air-sea exchanges beneath the 
tropical storm will be examined for a constant depth surface 
boundary layer interacting with the ocean surface. 

Several numerical simulations of hurricane development 
have been made with complicated multi-level models. These 
quire very short time steps, included some 
interaction with the ocean. The scope of the study is to 
model only the boundary layer beneath a specified vortex. 

The Riehl model is, therefore, used as a diagnostic model 
assuming that it aoueieeeat each time step. The steady-state 
vortex above the boundary layer is defined by Riehl's model, 
which is specified by air-sea interaction. 

Specification of the boundary layers under the circula- 
tion provides an opportunity to study the changes in the 
ocean as well as the vortex. An essential assumption carried 
over from Riehl was a steady, Pe rciitenll vortex which in 
early experiments, is stationary. Leslie and Smith (1970) 
point out that the importance of the surface boundary layer 


on the dynamics of a hurricane is due to the meridional 
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circulation of the vortex, providing a strong coupling be- 
tween the azimuthal and vertical velocity fields. However 
this study will use the simple boundary layer of Riehl. 
Several investigations have shown that areas of low sea- 
surface temperatures occur in the wake of hurricanes. Jordan 
(1964) analyzed sea-surface temperature for the months of 
typhoon seasons and reported marked cooling of the ocean 
surface with passage of the hurricane. Leipper (1967) 
analyzed oceanic temperature for hurricane Hilda of 1964. 
In the latter study, comparison of pre-storm merchant ship 
reports and post-storm data, gathered by a planned cruise of 
the area in the Gulf of Mexico where Hilda passed, indicated 
Significant cooling in the surface layer during the passage. 
Jordan suggested that cooling of the sea-surface temperature 
was due, in part, to upwelling of cooler water or vertical 
mixing caused by the winds. Jordan concluded that mixing 
was the more important factor. However, other authors such 
as Fisher (1958) have considered upwelling to be the most 
important. Leipper (1967) stated that in addition to up- 
welling of cooler waters and vertical mixing of the surface 
layer, horizontal advection of surface waters and the thermal 
transfer of energy from the ocean to the atmosphere are im- 
as factors in determining the observed oceanic thermal 
structures. Leipper concluded that the dominant process in 
the observed cooling is the Ekman suction resulting from 
wind mixing and that upwelling is the next most important. 


A theoretical two-layer ocean model developed by O'Brien and 


Jeu 





Reid (1967) to simulate circulation beneath a hurricane 
verified that upwelling was an important factor. 

In this research the oceanic model is an elementary 
attempt to simulate the ocean temperature field resulting 
from vertical mixing due to the air-sea exchange. Eventual 
refinements of the ocean model may include upwelling of 
cooler water and oceanic advection, and also movement of the 
hurricane over the ocean will allow prediction of the ocean- 
ic effects. Volgenau (1970) has incorporated the oceanic 
effect in terms of the hurricane heat potential (HHP) which 
1S deiined as the excess of heat content over that of water 
ae 26C. 

Those specific quantities examined in this study are the 
temperature and moisture fields of both the atmosphere and 
the ocean, the effects of the atmosphere upon the oceanic 
fields and vice-versa. Also the effect of varying initial 
maximum wind speeds will be examined. Eventual study will 
include the effects of a varying drag coefficient in the 
boundary layer. In general, this study is an attempt to 
create an interacting model of a hurricane boundary layer 
and a simple ocean model, which may be developed into more 
detailed descriptions of the effects of sea-air interaction 


ia hiaeari @anes. 
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II. DISCUSSION OF THE MODEL 


This study is an experiment with the steady-state, 
symmetrical hurricane as described by Riehl (1963); there- 
fore this discussion will follow relevant development in 
his paper. 

A symmetric, 17x17 rectangular grid is used in Cartesian 
coordinates with 15 nautical mile spacing. The center of 
the hurricane coincides with the center of the grid and the 
eye wall is linked to the first grid point outbound along a 
diagonal. Cartesian coordinates are used to link the sym- 
metrical hurricane to the model of an ocean within a rec- 
Weneusar refzions in the imetial experiments the hurri.cane 
is stationary, but eventually the hurricane will be moved 
with respect to the ocean model in a Cartesian coordinate 
system. Transformation from the arbitrary Cartesian grid 
to polar coordinates centered on the eye of the hurricane is 


part of the procedure. (See Fig. 2) 


A. REVIEW OF RIEHL'S MODEL 
1. Outflow Layer 
The initial parameters include the maximum tangen- 
tial wind component, Vp;, the radius of maximum wind, rj, 
along with the latitude of the storm. Riehl and Malkus 
(1961) have shown that conservation of absolute angular mo- 
mentum is a reasonable first approximation in the outflow 


layer. (See Fig. 1) This constraint leads to the following 
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expression for computing r,, the radius where v,=0, 


ie 3 \ % 1 
no= (2a ni /F ) (1) 


where f is the coriolis parameter. It is noted that by 
specifying the eye wall or maximum wind at a diagonal grid 
point, r, becomes a function of the grid spacing. Variables 
which are dependent on radius only are computed along a 
diagonal. Those variables include pressure, vertical motion, 
temperature and moisture advection, and large scale moisture 
PX - 

2. Inflow Layer 

Conservation of potential vorticity in the inflow 

layer was also demonstrated by Riehl and Malkus (1961). 
Matching the outflow and inflow layer at the eve wall leads 
Wo the relationship, 


y 
Uo It ~ = Const . (2) 


which is a good approximation to vg-r profiles found in 
mature hurricanes (see examples in Riehl, 1963). Inside 
the eye wall, solid body rotation is assumed to determine 
Vg: Thus given vg, and r, the tangential component of the 
wind can be determined from Eq. (2) at any radius. In par- 
ticular the tangential wind at the radius r, is used in the 
calculation of the initial radial component of the wind. 
Integration of the momentum equation over a shallow 
inflow layer of thickness ap where potential vorticity is 


assumed to be conserved, results in the following expression 









outflow layer = ©, 


Illustrating hurricane mass circulation 
Ve is tangential wind component, Geis 
equivalent potential temperature, r is 
radius 


Figure 1 





Illustrating polar coordinate system used 
ve is tangential wind component, v, is 
radial wind component, r is radius, Qis 
angular measurement 


Figume 2 
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for the radial wind component, 


Mn = ¥ /oP Ho poe | te. (*8 on Zin, (3) 


where g¢ is the acceleration of gravity and AP is the inflow 
thickness set at 100 mb. | 

The tangential component at the surface (Vg,) must 
also be specified to define a complete surface wind for 
calculation of the air-sea boundary fluxes. Therefore vg, 
was assumed to be 0.7 of the wind speed at the top of the 
boundary layer, which was specified by Eq. (2). The use of 
this coefficient is equivalent to specifying an effective 
drag coefficient based on the tangential component at the 
top of the boundary layer. Choosing a value of 0.7 reduces 
th 


ad 


)} by a factor of 0.49 for momentum 


@ 
CL 
Fs 
ap) 
0"3 
O 
© 
0) 
r+) 
Fy 
+4 
Q 
+4 
CD 
4) 
ct 
“oom 
ae 


exchange, and by 0.7 for the heat and moisture exchanges. 
As noted in the introduction the heat source depends line- 
arly on wind speed while momentum sink depends on the square 
of the wind speed. 
3. The Steady State Value of V 

Riehl proposed that a particular relation exists 
between the heat source and momentum sink in the interior of 
hurricanes, if a steady state velocity field is observed. 
This proposal leads to the following expression for the 


maximum wind speed, 


3 Man 


ae (issn) bE Ls ae i fie | (4) 


which satisfies both the tangential distribution from 
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momentum and from heat source considerations, where r, is 
the point where vg =25 m/sec. as chosen by Riehl. An ener- 


gy exchange potential parameter, E, is defined as, 


E = (Tw-Te lt L/Ce (Qw-Qa) (5) 


where Tw and Qw are the temperature (C) and the saturated 
specific humidity (gm/kgm) at the sea surface and Ta and Qa 
are the temperature and specific humidity of the atmosphere 
in the boundary layer. L represents the latent heat of 
vaporization (597.3 cal/gmw) and Cp is specific heat con- 


Seant for (.2k cal/gm deg). 


B, PREDICTION EQUATIONS 
i Atmespneric Prediction 

In the atmosphere the predicted variabies are the 
potential temperature and the specific humidity above the 
sea surface. These two parameters, along with pressure, 
determine the equivalent potential temperature Ge. Riehl 
(1963) has shown that Q@p is one factor in determining the 
maximum wind. 

The thermodynamic equation is used in the following 
forn, 


dot hr» att ds? 9) 


ae. ues + Wand 29 i 
where © is the potential temperature (K),UWis the vertical 
motion (mb/hr), P is pressure (mb) and Qs is sensible heat 


transfer. The first term on the right hand side of Eq. (6) 


is neglected due to symmetry. Potential temperature is 


Ly 





assumed to be constant in the boundary layer and thus could 
be created from the air temperature (Ta) field knowing the 
Beerece pressure. The effect of expansion cooling is in- 
cluded by using potential temperature in the thermodynamic 
equation. 
Moisture changes were determined in similar fashion 
using, | 
Qe. _ | 2m ay eB, , Hage Ou: 2 
»1t nr»re Addn >? 
where Qa is the specific humidity (gm/kgm) and Qe is the 
latent heat transfer. The first term on the right hand 
Side is neglected because of symmetry. 
Because the ocean model is in Cartesian coordinates, 
surface excnange valves are computed in rectangular coordi- 
nates. The transfer of sensible heat per unit volume is 


computed with the bulk transport equation in the forn, 


Qs /Ce 3 2 es Gime \ ut +ar™ (8) 


in units of cal/sq cm sec where Cp is specific heat, Ps is 
the atmospheric density above the sea surface, Tw is the 
temperature of the ocean, Ta is the atmospheric temperature, 
and Vite ay? represents the total Wind speed. Latent heat 
transfer or the rate of moisture input per unit volume is 
calculated by, 

@e/L = Cs (Qw-Qa) Vitra (9) 


in units of cal/sq cm sec where L is the latent heat of 
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vaporization, Qw is the specific humidity at the tempera- 
ture of the ocean surface (Tw) and Qa is the specific 
humidity of the atmosphere. The Bowen ratio is defined as 
the ratio of sensible heat to latent heat exchange and is 
an indication of the relative importance of these two 
mechanisms of energy transfer from the ocean to the atmos- 
phere. 

For use in Eq. (8) and (9), the wind is delineated 
into u and v Cartesian components. Cartesian components are 


transformed from the radial and tangential components of the 


wind as follows, 


Al = in Cos @ - Us 
(s- Mn SIN © + Ne 


Sin @ 


CoS © (10) 


J 


A 


ied 


counter-—clockwise direction 
from the x-axis so that u is positive for westerly winds and 
Varo positive for southerly winds. 

The equation of motion in polar coordinates, neg- 
lecting lateral friction, was used to determine the radial 


component of the wind, v , which was in balance with the new 


tangential wind component obtained from the following 


equation, 

diye , (Uo Un 1 tin =- 1 dP oes (iy) 

at H PARE ff OE 
In pressure coordinates, assuming a steady state, Eq. (11) 
yhebds, Oe 37 

a — AP 

f, Wo, As (12) 
dh ne 
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where 7 is the wind stress, )P/36@ is zero due symmetry of 
the storm, W je is zero because wis zero at the surface 
and We is 5 ar zero at the top of the boundary layer. 
A oe layer saeoe of 0.7 was Incorporated to account 
for the vertical profile, as was described in the inflow 


layer section, II A. The surface stress is defined to be 


proportional to the square of the surface wind, 


Tes = Ca fs We. / Cos & ) GiB) 


where o{ is the inflow angle measured relative to isobars. 
The stress 1s assumed to be zero at the top of the boundary 
layer. 

Computations of vertical motion and advection in 
Cartesian coordinates did not prove to be successful. This 
was because any values such as vertical motion or the ad- 
vection quantities, which involve x and y derivatives were 
not symmetrical with respect to the diagonals. Therefore, 
asymmetries resulted in the temperature and moisture fields 
and prevented long period predictions. Transforming to a 
polar coordinate system for the computation of vertical 
motion and the advection quantities conserved the symmetry 
_of the initial fields. Also, for asymmetric temperature 
and moisture fields, a small error is introduced by the 
neglection of the azimuthal advection. Due to the large 
tangential velocities wm a hurnecamess var hous properties 
tend to be symmetric with respect to the center and omitting 


tae Aagmmcnal advection introduces a mmen smaller error than 


the use of Cartesian. coordinates. 
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Vertical motion at the top of the boundary layer was 


computed from the continuity equation, in polar coordinates, 


de i de Ve " 3W ee) (iz) 
Lrwe rn IP 7 


where the change of the tangential wind with respect to 
angular change is zero in a symmetric storm. Vertical 
motion at the surface is assumed to be zero and the bounda- 
ry layer is chosen to be 100 mb thick as an approximation to 
the first kilometer of a tropical atmosphere. 
Ee 6(OcCeamm cc Prediction 

Predicted variables in the oceanic portion of the 
model are the sea-surface temperature (Tw) and the saturated 
specific humidity (Qw) at temperature Tw. A crude approxi- 
MeeioOn of the oceanic enerey loss centers on the concept of 
hurricane heat potential, HHP, defined by Volgenau (1970) 
as the excess of heat content over that of water at 26 C. 


A temperature profile as shown - pe = 20” 


loom mrs 29° 


@ 


22 
defines a hurricane heat potential of 3x104 cal/sq cm by, 


‘oe 
HHP = Ce (T- 26) x depTH , (15) 

where T is the mean temperature of the chosen ocean layer 

with a depth of 100 meters. The loss of energy by the 

ocean or the change in heat potential is equal to the sensi- 


ble and latent transfer of energy during each time step, 
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that is, 


AHHP = at x 3600.0 x(Ast+Qe) , (16) 


moore At is the prescribed time stepy After a time step 
the reduced magnitude of heat potential corresponds to a 

new mean temperature in the upper layer of the ocean and 
thus a sea-surface temperature can be calculated. This 
scheme accounts for sensible heating and evaporation but 
neglects radiation gains, horizontal divergence or upwelling. 
Since the saturation specific humidity field is computed 
from the Clapeyron equation, the evaporation from the sea 
surface is decreased as the surface temperature is de- 


creased. 


C. INTERACTION BETWEEN THE ATMOSPHERIC AND OCEAN MODELS 

In the atmospheric model the maximum wind speed deter- 
mined by the heat source is a function of the energy ex- 
change potential, E, as given by Eq. (4). A new mean value 
of E was calculated after each time step using the newly 
predicted atmospheric and oceanic parameters. However the 
ocean is effectively an infinite heat source and the result 
of this approach was to bring the atmosphere into Sc aee 
rium with the ocean, that is, E tended to zero. Consequent-— 
ly the maximum wind speed predicted from Eq. (4) also 
decreased to zero although very warm and moist air was 
available in the atmospheric boundary layer. 

In a second approach, the maximum wind speed was calcu- 


lated from the equivalent potential temperature at the eye 
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wall as follows, 


} 
Uer (m/sec) =14 4) ( Gei - ©.) F (ye) 


where Qe. is set at 350 K. The tangential wind component 
was computed at other grid points by the vg-r relationship 


as given by Eq. (2). 


D. INITIALIZATION OF SURFACE TEMPERATURES AND HUMIDITIES 
Integration of the hydrostatic equation for typical 


values of Ge ina hurricane produces a relation, 


P = -2.56 Oe a 


: J 


where P’ is the pressure departure from 1005 mb and 6 is 
the equivalent potential temperature departure from 350 K. 
This expression specifies eye wall pressure variation caused 
by the oceanic heat source. With the use of Eq. (17) and 
(18), the pressure and equivalent potential temperature at 
the eye wall may be specified for a given maximum wind 
speed. The pressure distribution along each diagonal was 
then computed from the eradient wind equation. The pressure 
at intermediate points throughout the symmetric field was 
determined by linear interpolation using pressure at the 
eye wall and the distribution of Ve with radius. 
Initialization of the temperature and moisture fields 
both in the boundary layer and at the ocean surface are. 
important parts of this model. In the initial experiments 
the sea-surface temperature, Tw, was initialized at a con- 


stant value of 30 C throughout the grid. The available 
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moisture at the ocean surface, Qw, was calculated from the 
integrated form of Clapeyron's equation for the saturation 


weeper pressure over water, 


hues 22s Tee ree (19) 
4.105 T 213 


where @€, is the saturation vapor pressure. Assuming the 





air in contact with the water is saturated, specific humidi- 


ty was calculated from, 


Ws = 622 Cs (20) 
te. Shes 


ia Uniitseof em/kem. 

Two different methods were used for initialization of 
the air temperature and moisture fields within the boundary 
jayer. One approach was to snecify the initial vaiues of 
the air temperature, Ta, and the specific humidity, Qa, in 
terms of the equivalent potential temperature at the eye 
wall which is related to the maximum tangential wind in 
Eq. (17). Specification of initial parameters at the eye 
wall from equivalent potential temperature is consistent 
with the scheme used in predicting the new maximum wind 
speed. This method involved the solution of two simulta- 
neous equations defining E and thé pseudo-equivalent tempera- 


ture. These equations were, 

@a= Qu - G/L fe = (w-Ta)] (21) 
amd, 

Gao = Ta SP Ly Tae 


L Ta (22) 


eh 





where Tae is the pseudo-equivalent temperature. Following 
Riehl, it was assumed that E may be treated as constant 

over the inner region of a hurricane and the equivalent 
potential temperature difference between the:‘eye wall and 
the outer regions of the storm is constant over the boundary 
layer depth Az. Thus the expression for the effect of the 
total heat source in increasing the equivalent potential 


temperature is, 


Gm: - Oc, = Cp E f wedt, (23) 


Az 


Effectively this expression implies bulk aerodynamic trans- 
fer of neat and moisture ina fixed boundary layer. The 
semeoabie and latent heat transferwwhere r<r, isv*erucial in 
determining the increased equivalent potential temperature 
required to generate hurricane force winds. The Ta of E 
was computed from Eq. (4). The adiabatic equivalent tempera- 


ture is related to the equivalent potential temperature by, 


Ge = Toe = (2h) 
? : 
By specifying the initial state from the equivalent poten- 
tial temperature a the eye wall, the effect of expansion 
cooling and sensible heat source at the surface are taken 
into account. Riehl suggested that Oe must rise above 
350-352K to allow deep convection which warms the column 
and lowers the central pressure to vaiues which can sustain 
hurricane force winds. 


Solution of Eqs. (23) and (24) gives an initial value 


of surface temperature and moisture at the eye wall. 


cae, 





Therefore expression Eq. (17) relating vg. to 6. applies 
only at the eye wall unless these values are specified 
elsewhere. To initialize the remainder of the Ta and Qa 
values it was assumed that the relative humidity varied 
linearly between the value determined at the eye wall anda 
value at the outer radius. The latter value was selected 
to be typical of undisturbed tropical air. For the remain- 
ing grid points simultaneous solution of the energy ex- 
change potential, E, equation and the Clapeyron equation 
including the specified relative humidity was used to 
determine the initial state. 

The Be seat experiment has shown that a second approach 
may be as effective as the previous solution of simultaneous 
equ@euzons, Py initializing the atmospheric temperature at 
an arbitrary value, say equal to the oceanic value, while 
the moisture was made consistent with the maximum wind speed 
through the equivalent potential temperature, it was found 
that the model came into balance with a steady-state hurri- 
cane after 10 hours. Comparison of the temperature and 
humidity fields produced by both approaches shown in Fig. 3 


and 4, demonstrates this favorable result. 
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Same as Fig. 3 except for specific humidity at the eye 
wall. 


Figure 4 


a | 





IIT, ATMOSPHERIC MODEL PROFILES 


Results of the basic model described in earlier sec- 
tions may be compared with the observed surface structure 
ef mature tropical storms. Specifications of the tangen- 
tial wind profile, vg , according to Eq. (2) has been shown 
by Riehl (1963) to be realistic, yet the remaining variables 
of the predicted hurricane need to be verified. Observa- 
tions of the surface characteristics of a hurricane have 
been summarized by Riehl (1954). For these comparison 
experiments the sea-surface temperature, Tw, was held fixed 
to model the steady-state atmospheric circulation. 

The effect cf the surface stress is to develop a cross- 
isobaric flow, the radial component of the wind (v,), shown 
in Figure 5. The radial wind component increases toward 
the eye wall, with a maximum value of v, which compares 
favorably with the velocity profile for a steady-state 
hurricane described by Riehl. Vertical motion is determined 
by the radial gradient of rv, in an axisymmetric storm. 

All values of W were negative with the strongest ascent 
concentrated near the center of the storm. Figure 6 illus- 
trates the area of strong upward vertical motion resulting 
from the inflow near the eyewall. Note that no strong radi- 
al inflow bands, which would represent the satellite- 
observed cloud formations spiralling toward the hurricane 


center, can exist in the axisymmetric vortex. 
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Same as Fig. 5 except for vertical motion in pressure 
coordinates. 


Figure 6 
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Although potential temperature of the surface air in- 
creases along the inflow trajectories, the surface tempera- 
ture decreases toward the storm center until a balance of 
expansion cooling and sensible heating is reached. The 
expected temperature pattern is verified by Fig. 7 with the 
maximum cooling of the surface air occurring at the eye 
wall. Specific humidity may be expected to increase along 
the inflow trajectories as the air in contact with the ocean 
acquires moisture through evaporation, which should increase 
with wind speed. Figure § shows a slight increase of mois- 
ture content predicted as the air moves inward. 

It has been observed that precipitation in hurricanes 
1s strongly concentrated near the inner core of a storm. 
Riehl (1954) notes the average rainfall rate at a fixed 
point will be on the order of 15-20 cm per day. For com- 
parison the large scale moisture flux due to vertical motion 
through the top of the boundary layer was converted to an 
equivalent depth of precipitation. The model predicts only 
3-4 cm per day, but Fig. 9 verifies the concentration of 
rainfall at the eye wall. 

The surface layer air acquires both latent and sensible 
heat as it moves inward toward lower pressure. These results 
Siewe chau tie transter of senstple heat from the ocean to 
the atmosphere is small compared to the latent heat transfer. 
In the initial state the air temperature, Ta, is approxi- 
mately equal to the sea-surface temperature, Tw, and the 


atmospheric specific humidity, Qa, is small compared ~o the 
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same aS Fig. 5 except for atmospheric temperature. 
Right hand scale represents the difference between 
oceanic and atmospheric temperature. 
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Same as Fig. 5 except for accumulated precipitation 
from large scale moisture flux after one hour. 
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Same as Fig. 5 except for sensible heat (Qs) and 
latent heat (Qe) exchange for the surface. Ocean 
temperature is fixed at 30C. 


Figure 10 
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saturation specific humidity at the ocean interface, There- 
fore latent heat transfer is larger than sensible heating. 
Figure 10 is a comparison of the heat transfer. The large 
addition of moisture due to evaporation may either be 
precipitated or stored in the column. The large precipita- 
tion at the eyewall is due to the fact that the moisture 
which is not precipitated along the radial inflow trajec- 


tory will be carried inward to the eye wall. 
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IV. MODIFICATION OF THE MODEL 


An encouraging aspect of the initial experiments was 
obtaining the expected atmospheric horizontal temperature 
profile. However, prediction of moisture content repre- 
sented a major shortcoming of the initial model. That is, 
the advection scheme plus the sence heat exchange pro- 
duced the expected temperature profile, but advection plus 
evaporation produced an excessive increase in moisture in 
the boundary layer. As shown in Fig. 12, surface relative 
humidity values at the eye wall exceeded saturated condi- 
tions after short time periods. Either the moisture trans- 
poms dwe to werticad metion in Ea. (7) was not removing 
enough moisture from the boundary iayer, or the evaporation 
from the sea was too large. 

A possible explanation may be that the surface exchange 
of moisture was overestimated. Robinson (1966) reports 
tCieeecic’ drag cocificient, C,, for the sensible and latent 
heat exchange at high wind speeds is less than the drag 
coefficient for momentum exchange by a factor of one-third. 
In the momentum exchange at the sea surface the pressure 
forces dominate over molecular processes at the interface. 
However near the air-sea interface molecular processes 
determine the rate of transport of heat and water vapor. 
Thus the transfer coefficient for momentum may differ from 


that for sensible heat and moisture exchanges which are 
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Figure 12 
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resumicted by molecular diffusion. To test this effect, 
computations were made with the sensible and latent heat 
exchange in Eqs. (8) and (9) reduced by a factor of one- 
third. Taking this into account served to decrease the 
surface air temperature and the moisture content as shown 
by Fig. 13. It should be noted, if the drag coefficient 
for momentum had also been reduced, the radial and vertical 
wind components would have been decreased. A reduced hori- 
zontal and vertical transport would have resulted which 
would have offset the favorable reduction in evaporation. 

Another explanation for the accumulation of moisture 
may be that the model includes no representation of the 
cumulus clouds or "hot towers" that characterize a hurri- 
cane eye wall. Such cumulus convection would transport 
additional moisture through the top of the boundary layer. 
Two approaches were used to reduce the excessive accumula- 
tion of moisture. The drag coefficient reduction suggested 
by Robinson was included in all experiments comparing these 
convective mechanisms. | 

Convective mechanisms were based upon relative humidity 
Erieerta, in each time step, a small amount (initially 0.5 
percent) of moisture was removed if the relative humidity 
exceeded a critical value. The amount of moisture removed 
could be varied with increasing relative humidity to achieve 
a balance between the moisture increase due to the hori- 
zontal advection and the convective moisture reduction. 
The amount of moisture removed by this scheme was accumu- 


lated and converted to an equivalent depth of precipitation. 


36 





30 


29 shilaa pinta: without Cd reduction 


with Cd reduction 


28 


27 


26 


25 





Wee HS. 


29 


28 eo 
without Cd a? 
reduction 





0 5 10 15 
TIME HRS, 


Evolution of atmospheric temperature (upper diagram) and 
specific humidity (lower diagram) with time for experi- 
ment with and without Robinson's C reduction. 


Figure 13 
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The first attempt to simulate cumulus convective mechanism 
is illustrated in Fig. 14. It was assumed that a convective 
transfer of moisture through the top of the boundary existed 
which was just sufficient to maintain the relative humidity 
at some critical value. 

In experiments to establish the critical value for use 
in this convective mechanism, it was found that the Qe at 
the eye wall was very sensitive to small variations. The 
percentage of specific humidity to be removed was varied 
along with the relative humidity criterion which initiates 
the mechanism. Figures 16 and 17 show that variation of 
either the percentage or relative humidity criterion caused 
a varying specific humidity content. The specific humidity 
1s crucial in determining the equivalent potential tempera- 
ture at the eye wall, Qe; , which specifies the maximum 
wind speed, Eq. (17). Small modification procedures pro- 
duced rather large variations in the intensity of the hurri- 
Cane. Therefore specification of the critical relative 
humidity essentially determined the steady-state hurricane. 

The first attempt at moisture reduction led to some 
improvement in the model. By preventing super-saturation, 
the moisture field was more representative. One of the 
objectives of the study was to achieve a steady-state condi- 
tion in the hurricane when the ocean temperature was held 
fixed. However it was found that regardless of the initial 
maximum tangential wind speed, the predicted hurricane had 


the same intensity, as shown in Fig. 18. For experiments 
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descent in the environment. 









Figure li 
upward Cu downward 
Vemtical weied ca | 
motion motion 


toss of Moisture 
from precipitation 


eye wall 


evaporation 


Convective moisture mechanism with Ta adjustment. 
Note that warming from the release of latent heat 
occurs in boundary layer. 
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Specific humidity with varying amounts of convective 
moisture reduction (fraction of 1% removal) for a 
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Specific humidity with varying relative humidity cri- 
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Maximum tangential wind speed (upper diagram) and 
specific humidity (lower diagram) for various initial 


Ve, values with a convective mechanism without a tempera- 


ture adjustment. | 
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with high initial wind speeds, and thus high specific hu- 
midity values near the eye wall, the removal of moisture by 
increased convection was not offset by an increased poten- 
tial temperature of the surface air. Therefore there was a 
Steady decrease in equivalent potential temperature at the 
eye wall, with a corresponding decrease in the wind speed, 
until a balance between the convective flux of moisture and 
the evaporation was achieved. This led to a second approach 
to offset the decreased specific humidity with an increase 
in the atmospheric temperature. 

Recall that the equivalent potential temperature, Qe, 
is conserved in a moist adiabatic process. In the second 
approach (see Fig. 15), the percentage moisture reduction 
is still utilized, but the process is assumed to be a psuedo- 
adiabatic process. This modification was accomplished by 
computing the equivalent potential temperature, Ge, from 
Eq. (24). The equivalent potential temperature is main- 
welimec aiber the reductionsot the specific humidity by 
increasing the potential temperature of the air to a value 
given approximately by, 

one -286 L Ga 
Tae - 2732 +0e/('S°) e Gt (25) 

where the previous value of Ta is used in the denominator. 
Yet Fig. 19 demonstrates that too much warming accompanied 
the moisture decrease. The equivalent potential temperature 
increased slowly, resulting in a corresponding growth in 


maximum tangential wind speed. Therefore a second 
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"steady-state hurricane" has been described fOr asdin ferent 
specification of the convective mechanism. 

These experiments to reduce the predicted excess mois- 
ture point out the sensitivity of the model. Evidently a 
combination of moisture transport out of the boundary layer, 
along with some precipitation warming within the layer, is 
necessary to produce a steady-state hurricane. Experimen- 
tation 1s necessary to relate the combination of these two 


effects to the intensity of the storm. 
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V. OCEAN MODEL PROFILES 


Two nearly steady-state models of the atmospheric flow, 
presented in the modification of the model section, allow 
comparison of the effect of hurricane circulation upon the 
ocean profile. The ocean temperature was allowed to vary, 
due to wind mixing. 

A comparison of hurricane heat potential (HHP) and the 
changes in sea-surface temperature, Tw, will be made to | 
show the effect of the two hurricanes. Hurricane heat 
potential is a function of the sensible and latent heat 
transier..as given byeEqs. (8) and (9). The change in this 
parameter represents the total energy loss experienced by 
the ocean due to the air-sea boundary fluxes. The greatest 
loss of heat potential occurred near the eye wall where the 
latent and sensible heat transfer-was the largest, as shown 
by Figs. 20 and 22. The accumulated energy loss increased 
with time with a loss of 1600-1900 cal/sq cm per day being 
accounted for by this approach. This energy loss is small 
compared to an atmospheric heat balance calculation of 
Malkus (1962), in which the total* énergy loss approached 
3000 cal/sq cm per day. In the present model the hurricane 
heat potential calculation does not include radiation, or 
advection and upwelling effects. 

Comparison of Figs. 20 and 22 indicates that the heat 


potential loss near the eye wall was less for the atmospheric 
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Hurricane heat potential loss for varying daily 
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Sea-surface temperature loss for varying daily 
observations for a convective mechanism without a 
temperature adjustment. 
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model with a convective mechanism with the temperature ad- 
justment. This result was some what surprising since the 
wind speeds were larger in this model. The sensible and 
jatent heat transfers near the eye wall were smaller be- 
cause the air-sea temperature and moisture differences were 
smaller. In this region where convective heating was im- 
portant, the temperature tends to increase through precipi- 
tation warming. As the sea-surface temperature, Tw, 
decreased and the air temperature, Ta, increased, the 
sensible heat transfer was smaller. A similar explanation 
can be given for a reduction in evaporation. That is, the 
specific humidity, Qa, was larger (for the same critical 
relative humidity) in the adjusted temperature case because 
the caturation limit increases with temperature. Away from 
the eye wall, the wind speed became the dominant factor in 
the sensible and latent heat transfer equations. As shown 
in Fig. 22, the change in HHP was. larger at these radii for 
the convective mechanism with temperature adjustment. 

The simple model of calculating the sea-surface tempera- 
ture from the loss of hurricane heat potential (HHP), 
Eqs. (15) and (16), resulted in a temperature field which 
cooled with time. Since the coolting occurs over a deep 
layer a nearly uniform cooling rate of .3 — .4C0 per day was 
noted near the eye wall region. This cooling rate is small 
compared with observed values of 1 - 2C cooling, as might be 
expected since radiation, horizontal advection in the ocean 


and upwelling effects have not been included. The resulting 
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sea-surface temperature fields point out the same features 
when comparing the two different convective mechanisms as 
were noted in the hurricane heat potential discussion. 
Evolution of the maximum tangential wind speed over an 
extended period of integration is shown in Fig. 24. As 
noted previously the maximum tangential wind speed had a 
lower value without a temperature adjustment than when the 
surface temperature was adjusted for the effect of precipi- 
tation heating within the boundary layer. When the ocean 
temperature was allowed to vary, the resulting maximum wind 
speeds were slightly lower even in the first few hours, 
Pe. 2h. As ohne ocean surface cooled, the sensible heat 
transfer decreased, and produced slightly cooler atmos-— 
pheric temperatures at the eye wall. Again, restrictions 
on the relative humidity value reduced the specific humidity, 


resulting in a lower tangential wind component. 
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Comparison of maximum tangential wind speed with 
sea-surface temperature fixed (Tw fixed) and sea- 
surface temperature allowed to vary (Tw free) for 
the model incorporating a convective mechanism with 
and without temperature adjustment. 
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VI. CONCLUSION 


This research utilized the steady-state hurricane model 
proposed by Riehl in experiments with a simple interacting 
ocean model. The purpose was to examine the effect of the 
air-sea exchange ina fixed boundary layer on the hurricane 
vortex and on the sea-surface temperature. 

The model predicted fields of hurricane parameters which 
compare with observed surface features of mature tropical 
storms. The radial inflow profile within the boundary layer 
ws realistic, with a corresponding vertical motion field 


Similar to observed values. Because the model does not have 


ah 


detailed boundary laver, the solution of the radial wind 
component due to surface stress is not very sophisticated. 
More resolution in the boundary layer and inclusion of 
differences in surface roughness would possibly produce a 
better radial wind component solution. 

The simple representation of temperature advection plus 
sensible heat exchange produced the observed inward decrease 
in temperature in the atmospheric boundary layer. However 
an excessive accumulation of moisture occurred in the bounda- 
ry layer. Two attempts to correct for this situation have 
shown that small changes in the surface moisture and tem-— 
perature values strongly affect the maximum wind speed 
predicted. The proportion of the water vapor transported 


out of the top of the boundary layer without adding latent 
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heat to the boundary layer appears to be an important factor 
in determining the intensity of a steady-state hurricane. 
Although the ocean model is quite simple, realistic 
horizontal oceanic temperature profiles with respect to the 
hurricane center were obtained. However the magnitude of 
changes in oceanic heat content and the sea-~surface tempera- 
tures were smaller anew expected in actual hurricanes, since 
advection and upwelling effects were not included in the 
ocean model. A steady-state cannot be reached in a station-— 
ary storm because of the decreasing sea~surface temperature 
field. Provision for movement of the hurricane with respect 
to the ocean would allow prediction of the intensity of the 
atmospheric system and the oceanic effects. Improvement of 
the oceanic model, to include a combination of vertical 
mixing and upwelling, would also give @ more realistic view 


of the consequences of air-sea interaction below a hurricane. 
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